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Abstract 

We have performed a systematic study of the electronic structures of BiMeOs {Me = Sc, Cr, 
Mn, Fc, Co, Ni) series by soft X-ray emission (XES) and absorption (XAS) spectroscopy. The 
band gap values were estimated for all compounds in the series. For BiFcOs a band gap of --^0.9 
eV was obtained from the alignment of the O Kq XES and O Is XAS. The O Is XAS spectrum of 
BiNiOs indicates that the formation of holes is due to a Ni^+ valency rather than a Ni^+ valency. 
We have found that the O Kq, XES and O Is XAS of BiMeOs probing partially occupied and 
vacant O 2p states, respectively, are in agreement with the O 2p densities of states obtained from 
spin-polarizcd band structure calculations. The O XES spectra show the same degree of Bi 
6s — O 2p hybridization for all compounds in the series. We argue herein that the stereochemical 
activity of Bi 6s lone pairs must be supplemented with inversion symmetry breaking to allow 
electric polarization. For BiMnOa and BiFeOa, two cases of multiferroic materials in this series, 
the former breaks the inversion symmetry due to the antiferromagnetic order induced by particular 
orbital ordering in the highly distorted perovskite structure and the latter has rhombohedral crystal 
structure without inversion symmetry. 
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I. INTRODUCTION 



Multiferroics, discovered almost 50 years ago^, are materials that simultaneously pos- 
sess two or three degrees of freedom: (anti)ferromagnetism, (anti)ferroelectricity, and/or 
ferroelasticity, which allow both charge and spin to be manipulated by applied electrical 
and magnetic fields^'^. These materials are promising for various technological applications 
such as information storage, spintronics, and sensors. There are many different classes of 
multiferroic systems known today, for instance: the i^MnOs family {R = Dy, Tb, Ho, Y, 
Lu, etc.), the i^MngOs family {R = Nd, Sm, Dy, Tb), and the BiMeOg family (Me = 
Mn, Fe). These materials have complex structures with many atoms per formula unit, and 
more than one formula unit per unit cell. The large number of interatomic interactions 
makes distinguishing the mechanisms responsible for multiferroicity a challenging task. The 
origin of these phenomena and the nature of the coupling between the magnetic, electric, 
and structural order parameters are not well understood. Despite many electronic struc- 
ture calculations (see review article [sl and references therein) only few experimental spectra 
for selected compounds (YMnOs and BiFeOs) have so far been obtained^*^. In order to 
get deeper insight into the nature of multiferroicity we have performed a systematic study 
of electronic structure of the perovskite-like multiferroics (BiFeOs and BiMnOs) and the 
BiMeOa (Me = Sc, Cr, Co, Ni) related compounds using synchrotron excited soft X-ray 
emission and absorption spectra. The experimental results are compared with specially 
performed electronic structure calculations. 

The paper is organized as follows: the details of sample preparation and X-ray measure- 
ments are presented in Section II. The crystal structure of different compounds and their 
basic properties are summarized in Section III. Results of X-ray measurements for the whole 
series of BiMeOs {Me = Sc, Cr, Fe, Co, Ni) compounds and comparison with electronic 
structure calculations are collected in Section IV. 

II. EXPERIMENTAL DETAILS 

All the samples were synthesized using a high-pressure high-temperature method. Start- 
ing mixtures were placed in Au capsules and treated at 6 GPa in a belt-type high-pressure 
apparatus at different temperatures (heating rate was about 140 K/min). After heat treat- 
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ment, the samples were quenched to room temperature (RT), and the pressure was slowly 
released. BiCrOa was prepared from mixtures of Bi203 (99.99 %) and Cr203 (99.9 %) with 
an amount-of-substance ratio of 1:1 at 1653 K for 60-70 min^. BiMnOs was prepared from 
mixtures of Bi203 and Mn203 with an amount-of-substance ratio of 1:1 at 1383 K for 60-70 
min^. Single-phased Mn203 was prepared by heating commercial Mn02 (99.99 %) in air at 
923 K for 24 h. For the preparation of BiScOs, stoichiometric mixtures of Bi203 and SC2O3 
(99.9 %) were dried at 873 K for 8 h and then treated at 1413 K (at 6 GPa) for 40 miiM. 

For the preparation of BiFe03, stoichiometric mixtures of Bi203 and Fe203 (99.9 %) 
were first annealed at ambient pressure at 1073 K for 2 h. This procedure gave a mixture of 
BiFe03 (about 70 weight %), Bi25Fe039, and Bi2Fe409. The resulting mixture was treated 
at 1273 K (at 6 GPa) for 1 h. After the high-pressure treatment, single-phased BiFe03 was 
obtained. BiCo03 was synthesized from stoichiometric mixtures of Bi203, C03O4 (99.9 %), 
and KCIO3 at 1483 K (at 6 GPa) for 1 h, and BiNi03 from Bi203, NiO (99.9 %), and KCIO3 
at 1483 K (at 6 GPa) for 1 h. The resulting samples were reground, washed with water, 
dried, and re-pressed at about 1 GPa at room temperature. 

The X-ray emission spectra (XES) were measured at beamline 8.0.1 at the Advanced 
Light Source (ALS) at the Lawrence Berkeley National Laboratory^^. The X-ray absorption 
spectra (XAS) were measured at the Spherical Grating Monochromator beamline at the 
Canadian Light Source (CLS) at the University of Saskatchewan^. The O Is XAS spectra 
were measured in the total fluorescence yield (TFY) mode, which provides more bulk sensi- 
tivity than electron yield methods do. The O XES was excited near the O Is ionization 
threshold to suppress the high-energy satellite structure. The instrumental resolving power 
(E/AE) was approximately 10^ for the XES measurements and approximately 5x10'^ for 
the XAS measurements. 

III. CRYSTAL STRUCTURE AND PHYSICAL PROPERTIES OF BIMEO3 
COMPOUNDS 

Of all the multiferroic compounds, the most extensive study has been devoted to the 
perovskite-type and related oxides. Below we give a short summary of crystal structure and 
basic physical properties of the BiMe03 {Me = Sc, Cr, Mn, Fe, Co, Ni) compounds. The 
space groups and lattice constants for the whole series are collected in Table [B The selected 
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bond lengths are shown in Tables [TTl 11111 and ITVl 

The structure of BiScOs has monoclinic symmetry with the space group C2/c^^. BiScOa 
is nonmagnetic and most of the interest in this compound is in the ferroelectric properties of 
the BiScOa-PbTiOa system. For the well known piezoelectric material Pb(Zr,Ti)03 (PZT) 
the Curie temperature T^ at the morthotropic phase boundary between the rhombohedral 
and tetragonal ferroelectric state is 386° C. It has been shown that for BiScOs-PbTiOs the 
T(7 reaches 450° while the piezoelectric coefficients in bulk are comparable to those of 
commercial PZT. 

BiCrOs crystallizes in an orthorhombic structure above 420 K with the space group Prima 
and lattice parameters a = 5.54568(12) A, b = 7.7577(2) A, c = 5.42862(12) A at 490 K^. 
A structural phase transition from an orthorhombic to a monoclinic structure occurs at 420 
j^_8ji4^ii6 Between 420 K and 7 K BiCrOs has a monoclinic structure with the space group 
C2/c and a = 9.4641(4) A, b = 5.4790(2) A, c = 9.5850(4) A, p = 108.568(3)° at 7 K^. A 
long-range G-type antiferromagnetic order with weak ferromagnetism develops below Tn = 
109 K and does not change down to 7 K®. Four anomalies of magnetic origin were observed 
at 40, 75, 109 and 111 K^-. Magnetic moments of Cr^+ ions were found to align along the 
monoclinic b axis in a similar manner to the direction of magnetic moments of Mn^"*" in 
BiMnOa^'''^^. The magnetic structure of BiCrOs was first predicted from ab initio electronic 
structure calculations^^. 



TABLE I: Crystal structure and lattice parameters for BiMeOs compounds. 





BiScOg 


BiCrOs 


BiMnOs 


BiFeOs 


BiCoOs 


BiNiOs 


space group 


C2/c 


C2/c 


C2/c 


R3c 


P4mm 


PI 


a, [A] 


9.8899(5) 


9.4641(4) 


9.5415(2) 


5.58102(4) 


3.71990(7) 


5.3852(2) 


b, [A] 


5.8221(3) 


5.4790(2) 


5.61263(8) 


5.58102(4) 


3.71990(7) 


5.6498(2) 


c, [A] 


10.0462(5) 


9.5850(4) 


9.8632(2) 


13.8757(2) 


4.71965(15) 


7.7078(3) 


a, [°] 


90 


90 


90 


90 


90 


91.9529(10) 


A n 


108.300(3) 


108.568(3) 


110.6584(12) 


90 


90 


89.8097(9) 


7, [1 


90 


90 


90 


120 


90 


91.5411(9) 



The bismuth manganite (BiMnOs) has a highly distorted perovskite structure and has 
been regarded as a multiferroic material. The ferroelectricity has been analyzed within 
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TABLE II: Selected bond lengths in BiScOs, BiCrOs and BiMnOs, in A. A = ^ Ef (^^)^ and 
lav = J2i where k is the i-th bond. 







J3iurU3 


oiiVlnU3 


Mei-UZxz 


Z.UoO 


i.yooUo 


l.yUooo 


iwel-Uixz 


z.iuyyu 


i.yyioo 


Z.199UO 


Mel-Udxz 


2.15779 


1.9781)5 


1.98549 


L\[iviei-(j ) 


z.UX lU 


U.Uo X iU 


/ .z X lU 


iWez-Uoxz 


z.uybu { 


1.9 ^o9U 


I.94I0I 


iWez-Uixz 


z.iiooo 


1.99ZDO 


1.9z4Ul 


iMfez-Uzxz 


1 '7/1 


z.Ulo9o 


z.z41 /^4 


A(Mez-U ) 


U.bx iU 


U.ox ID 


ol.oX lU 


Bi-U3 


1 r r" 

z.loooz 


2.25097 


2.24562 


Bi-(Jz 


z.iyziz 


2.23752 


2.21778 


Bi-01 


2.24633 


2.32584 


2.23928 


Bi-01 


2.55437 


2.44659 


2.46625 


Bi-03 


2.89682 


2.65246 


2.70996 


Bi-02 


3.01816 


2.78741 


2.83703 


A(Bi-O) 


18.7x10^3 


7.0x10^2 


9.9x10^3 



TABLE III: Selected bond lengths in BiFeOa and BiCoOs, in A. 



Fe-0x3 


1.97139 


Co-01 


1.71781 


Fe-0x3 


2.08310 


Co-02x4 


2.01593 


A(Me-O) 


8x10^^ 




37x10"^ 


Bi-0x3 


2.30512 


Bi-02x4 


2.25155 


Bi-0x3 


2.51513 


Bi-01x4 


2.79852 


A(Bi-O) 


1.9x10-3 




11.7x10-3 



first-principles electronic structure calculations^", and attributed to the chemical activity 
of the Bi 6s^ lone pairs^^. BiMnOs is the only ferromagnet among the discussed BiMeOs 
compounds with a Curie temperature above 100 K. The largest saturation magnetization 
has been reported to be 3.92 /i^ per formula unit^, which is close to the 4 /i^ expected for the 
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TABLE IV: Selected bond lengths in BiNiOs, in A. 



Nil -09-5^9 


1 Q89Q4 


Ni9-r)3v9 


9 n4Q5Q 


Ni3-05v9 


1 Q8n75 


Ni4-01'ir9 


9 05408 


Nil -0^x9 




Ni2-04x2 


9 n7nQ7 


Ni3-02x9 


9 in6fi9 


Ni4-Ofix2 


2 10Q27 

Z; . J. VJl/Zj ( 


Ni1-04x9 


9 98fiQ7 




2.13415 


Ni3-Ofix9 


2.12071 




9 907fiQ 


A(Mp-C)'\ 


46x10^^ 




3x1 n^^ 




Qx 1 0^'^ 




Qx10~^ 


Ril D"^ 


9 91 ^1 1 


Ri9 ni 












Ril 


9 9fi^zL1 


Ri9 04 












Ril 09 


9 "^9401 


Ri9 OR 












Bil-03 


2.41745 


Bi2-02 


2.13787 










Bil-01 


2.52033 


Bi2-05 


2.24213 










Bil-04 


2.67658 


Bi2-06 


2.25756 










A(Bi-O) 


4.3x10^3 




1.7x10^3 











ferromagnetic state. The ferroelectric hysteresis loop has been observed in polycrystalline 
and thin film samples of BiMnO^^, although the measured ferroelectric polarization was 
much smaller (about 0.043 /iC/cm^ at 200 K) than the one obtained in the first principle 
calculations for the experimental noncentrosjnumetric structure (about 0.52 yuC/cm^)^^. 

BiMnOs undergoes two phase transitions at temperatures 474K (monoclinic-monoclinic) 
and 770K (monoclinic-orthorhombic}^i^i2^. According to early experimental data BiMnOs 
was considered to have noncentrosymmetric C2 space group below 770 K—'^. Recently the 
crystal structure of BiMnOa was re-examined by Belik et al.^ and confirmed by neutron 
powder diffraction experiments by Montanari et alM. The new experiments reveal that 
BiMnOs below 770 K has a centrosymmetric C2/c space group with parameters given in 
Table [H A structural optimization performed using modern methods of electronic structure 
calculations has shown that the noncentrosymmetric C2 structure, which had been reported 
earlier converges to the new total energy minimum corresponding to the C2/c structure with 
zero net polarizatiort^'^. 

Since the C2/c structure of BiMnOs has inversion symmetry the hypothesis that electric 
polarization arises due to bismuth lone pairs is no longer valid. The magnetic mechanism 
of inversion symmetry breaking was considered recently in work 29. It was shown that the 
peculiar orbital ordering realized below 474 K gives rise to ferromagnetic (FM) interactions 
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between nearest-neighbor spins which compete with longer-range antiferromagnetic (AFM) 
interactions. The solution of the low-energy model for 3d states in BiMnOs revealed that 
the true symmetry is expected to be Cc. The solution of the realistic model indicates a non- 
collinear magnetic ground state, where the ferromagnetic order along one crystallographic 
axis coexists with the the hidden AFM order and a related ferroelectric polarization along 
two other axes^S. 

The perovskite BiFeOs is ferroelectric with Tc = 1103 K and antiferromagnetic with T^v 
= 643 K and a canted spin structure^-. The bulk single crystal has a rhombohedrally dis- 
torted perovskite structure with the space group R3c^ and lattice parameters presented in 
Table [H The G-type collinear antiferromagnetic spin configuration has been modified by 
subjecting it to a long-range (620 A) modulation leading to a spiral modulated spin struc- 
ture^. The spontaneous polarization of a single crystal is 3.5 /iC/cm^ along (001) direction 
and 6.1 /xC/cm^ in (111) direction at 77 K— . This value is significantly smaller than spon- 
taneous polarization of lead titanate (80 to 100 ;uC/cm^ with the Tc ~763 K). However the 
heteroepitaxially constrained thin films of BiFeOa display a room-temperature spontaneous 
polarization of 50-60 yuC/cm^, an order of magnitude higher than in the bulk^^. The spin 
polarized first-principles calculation within local spin density approximation (LSDA) using 
the pseudopotential VASP package with the optimizated lattice parameters for the bulk 
rhombohedral phase results in polarization of 6.61 yuC/cm^— in excellent agreement with 
experiment. The thin films were shown to have tetragonal-like structure. For the P^mm 
symmetry and lattice parameters measured for thin film the Barry-phase calculation yields 
a spontaneous polarization of 63.2 yuC/cm^-~, consistent with experimental data, but it was 
revealed that small changes of lattice parameters can lead to a dramatically different polar- 
ization. The magnetoelectric coefficient {dE/dH, E - electric field, H - magnetic field) was 
measured to be 3 V/cm-Oe at zero field^^. Later on the values of remnant polarization were 
increased to 55 /xC/cm^ for (001) films, 80 /iC/cm^ for (101) films and about 100 /iC/cm^ for 
(111) films^. The BiFeOs films grown on (111) have the rhombohedral structure as single 
crystals, whereas films grown on (101) or (001) are monoclinically distorted^^. The high- 
est remnant polarization ever measured for a ferroelectric material, 146 /iC/cm^, has been 
reported for BiFeOs thin films with tetragonal crystal structure in work^"^. A wide range 
of measured polarization values were shown to be in consistent with the modern theory of 
polarization^i^i^, which recognizes that polarization is a lattice of values rather than a 
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single vector—. 

BiCoOs is isotypic with BaTiOs and PbTiOs and has a tetragonal crystal structure with 
the lattice parameters a = 3.71990(7) A, b = 4.71965(15) A, c/a = 1.269 at 5K^^ BiCoOg 
is an insulator with T^v = 470 K. It has C-type antiferromagnetic order where the magnetic 
moments of Co^+ ions aligning antiferromagnetically in the ab plane and antiferromagnetic 
jab layers stack ferromagnetically along the c axis'*^. The high-spin configuration of Co^^ 
ions (S = 2) has been reported^. The magnetic moments are m = 3.24(2)yUB for T = 5K 
and m = 2.93{2)^b at room temperature^. A reduction of the observed magnetic moment 
compared to the expected value 4/xb may be ascribed to the covalency of Co-0 bonds. The 



42 were 3.14 for Bi and -2.13 for 02 corre- 



bond valence sums (BVS) at 300K obtained in 
sponding to the oxidation states -|-3 and -2, respectively. The BVS values were 2.68 for Co 
and -1.57 for 01 indicating the covalency effects. The spontaneous polarization in BiCoOs 
of 179 ;uC/cm^ has been predicted from first-principle Berry-phase calculations^. The ex- 
perimental observation of a ferroelectric hysteresis loop is problematic since the resistivity 
appears too low for the large applied electric field. Therefore it was proposed that BiCoOs 
should be regarded as a pyroelectric rather than a ferroelectric material (in pyroelectrics 
the dipole moments can not be reoriented by external electric field)^^. The C-type anti- 
ferromagnetic order with a reduced magnetic moment of 2.41 fiB have been also predicted 
from first-principle calculations^. Spin-polarized calculations with C-type antiferromagnetic 
order for BiCoOa predicts an insulating ground state with an energy gap of 0.6 eV. 

BiNiOa has been found to be an insulating antiferromagnet {Tn = 300 K)^ with a heavily 
distorted triclinic symmetry PI. The lattice constants are shown in Table [H X-ray powder 
diffraction data has revealed that Bi ions were disproportionately weighted towards Bi^^ 
and Bi^+ and therefore the oxidation state of the Ni ion was +2 rather than the expected 
_l_244j45 rpj^^ Curie constant estimated from the magnetic susceptibility of BiNiOa is close to 
that expected for S = 1 rather than for S = 1/2 system. This fact as well as BVS confirms 
the divalent nature of Ni. The electronic structure of BiNiOs has been performed by the 
full-potential method within LDA-I-U { U = 8 eV, J = 0.95 eV) approximation with the 
G-type antiferromagnetic spin configuration^^. At ambient pressure a insulating solution 
with the charge-transfer gap of 0.8 eV was obtained in agreement with the value of 0.675 eV 
estimated from the electrical resistivity^. From the powder neutron diffraction study it was 
found that the valence state of BiNiOa changes under pressure^. Both neutron diffraction 
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measurements and BVS show that under pressure the charge disproportionate state melts 
leading to the simultaneous charge transfer from Ni to Bi, so that the high-pressure phase is 
metallic Bi^+Ni^'^03. This transition takes place at 4 GPa pressure and structure changes 
to the GdFeOa-type with the Pbnm symmetry^^. First-principle calculations also reproduce 
the metallic character of high-pressure phase. 

IV. RESULTS AND DISCUSSION 
A. Results of LSDA calculation 

The electronic structure of BiMeOs (Me = Sc, Cr, Mn, Fe, Co, Ni) series was calcu- 
lated within the local spin density approximation (LSDA)^i^ with a linearized muffin-tin 
orbitals basis (LMTO)^ using the STUTTGART TB-LMTO program (version LMT047). 
For all compounds the experimental atomic positions and lattice constants shown in Table [T] 
were used. For the spin-polarized calculations we used the experimentally determined mag- 
netic structures described in previous paragraphs (for BiFeOs the G-type AFM order was 
assumed). The resulting energy gaps (E^"''^') as well as magnetic moments are shown in Ta- 
ble |V] along with the available experimental data. For BiScOs, BiCrOs, BiFeOa and BiCoOa 
LSDA gives an insulating solution. For BiMnOa LSDA results in a half-metallic solution 
although the material is known to be an insulator. The metallic state of BiNiOa obtained 
with LSDA calculations is also in contrast with experimental observations. The failure of 
LSDA is due to the improper treatment of Coulomb correlations. It is well known that 
LSDA-I-U improves the description of correlation effects in transition metal oxides due to 
treating the strong Coulomb repulsion between locahzed d states by adding a Hubbard-like 
term to the effective potential^. This calculation requires the Hubbard parameter U and 
exchange interaction J. Although these parameters can be calculated by constrained LDA 
method?^ at present work we just take typical values of f/ = 3 eV, J = 1 eV for BiMnOa 
and f/ = 8 eV, J = 1 eV for BiNiOs. 

The total and partial densities of states (DOS) obtained for BiMeOs series within LSDA 
and LSDA-I-U (for BiNiOs and BiMnOa) are presented in Figures [T] and O One can see 
that the valence band for all BiMeOa compounds is formed by the Me 3d states hybridized 
with the O 2p states. The low lying states at energy about -10 eV comes from Bi 6s states, 
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the so-called "lone pair" . These states are only slightly hybridized with 2p states of oxygen. 

TABLE V: The calculated (E^"'^ ) and estimated from XES and XAS spectra (A^) values of 
energy gap in comparison with the experimental data taken from literature (E^** ). The calculated 
magnetic moments on the Me ion (uiMe) along with experimental estimations (m^|* ) for BiMeOa 
compounds. The different values of magnetic moments in case of Me = Cr, Mn, Ni correspond to 
the nonequivalent Me atoms in the unit cell. 



Compound 


^Calc. [gY] 


A, [eV] 


E^^*- [eV] 


mMe [fJ-s] 




Config. Me^+ 


BiSc03 


3.3 


2.6 









(f 


BiCrOs 


0.88 


1.4 




2.63, 2.65 


2.558 




BiMnOs 


0.33 


0.9 


insulator^^ 


3.65, 3.64 


3.2I8 




BiFeOs 


0.51 


0.9 


1.3*^-, 2.5^ 


3.54 


3.7532 




BiCoOg 


0.72 


1.7 


insulator^^ 


2.41 


3.24^ 


d^ 


BiNiOa 


1.23 


1.1 


0.675^4 


1.7; 1.67 


1.76^5 

















B.CrO. ■ 



















Energy, eV 



^ , 1 , , , , 1 , , , , 
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[ Ml/ 


i JWJl'l\: , , ,- 


- — 


'uiCoO,- 


















-10 -5 5 
Energy, ev 



FIG. 1: LSDA (LSDA+U with [/ = 3 eV, J = 1 eV for BiMnOa and [/ = 8 eV, J = 1 eV for 
BiNiOa) partial Me 3d and O 2p density of states for BiMeOs {Me = Sc, Cr, Mn, Fe, Co, Ni) 
compounds. The Fermi level corresponds to zero. (Colour in on-line version.) 



BiScOs, BiCrOs, and BiMnOs all have a monoclinic structure and can be compared 
directly. The bondlengths for these three compounds are summarized in Table [Tll Although 
the space groups of these compounds are the same, the MeOe octahedra distortions are 
different for each and most pronounced in the case of BiMnOs. The parameter A in Table HTl 
indicates the degree of MeOg octahedra distortion. The largest A corresponds to BiMnOs 
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while in BiCrOs the CrOe octahedron is almost undistorted. The shortest Bi — O bond 
lengths also shrink from Sc to Mn. The oxygen atoms surrounding Bi do not form regular 
polyhedra so we take the six nearest neighbors and estimate the variations between the 
compounds. All three compounds display rather strong deviation from average Bi — O bond 
but the largest one is found for the BiScOs system. BiFeOa and BiCoOa have different 
structures. Five oxygen atoms surrounding the Co ion form a pyramid (see Table IIIip . 
This pyramid is contracted so the distortion parameter is rather large. In BiFeOa the FeOe 
octahedron is less distorted than in BiMnOs system. The oxygen atoms are distributed 
around Bi in BiFeOa almost evenly while in BiCoOa rather irregular arrangement of oxygen 
atoms can be seen from Bi — O bondlengths shown in Table IIIII In Table IIVI the distances 
between Bi — O and Ni — O for BiNiOs are shown. The low symmetry of BiNiOa leads to 
the four nonequivalent Ni atoms and two nonequivalent Bi atoms in the unit cell. Among 
four Ni octahedra formed by oxygen atoms the one around Nil is the most distorted. The 
oxygen neighborhood around the Bi atoms is rather uniform. 

For BiScOs and BiCrOs an insulating solution was obtained with LSDA calculations (see 
Figure [1]). BiScOs has a formal electronic configuration and calculated energy gap of 3.3 
eV. The valence band is formed by O 2p states while the conduction band is composed of 
Sc 3d states hybridized with O 2p states. BiCrOa has three electrons occupying t2g orbitals 
(see Figured]). The energy gap of 0.88 eV and magnetic moment of 2.63 and 2.65 fiB (for 
two nonequivalent Cr atoms) was obtained in spin-polarized calculation in good agreement 
with the value 2.55 fiB obtained in neutron diffraction experiments^. Both the valence and 
conduction bands are formed mainly by Cr 3d states with small admixture of O 2p states. 
LSDA+U produces the insulating state for ferromagnetically ordered BiMnOa in agreement 
with experiment (see Figured]). The calculated magnetic moments for two nonequivalent 
Mn ions are 3.64 and 3.65 fiB which is in agreement with experimental values shown in 
Table |Vl Both the O 2p and the Mn 3d bandwidths are wider than those of BiCrOs but 
the main contribution to the valence band comes from Mn 3d states. There is also a strong 
overlap of O 2p and Mn 3d bands. 

LSDA calculations result in an insulating solution for BiFeOa with the energy gap of 0.51 
eV. This is much smaller than experimental estimates^*^. However the magnetic moment 
value of 3.54 is in good agreement with experiment (see Table |V]). The calculated DOS 
agrees well with the results of previous ab initio calculations^. Within LSDA the high spin 
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(HS) state for the configuration of Fe ion was obtained as shown in Figure [H Both the 
valence and conduction bands are formed predominantly by Fe 3d states hybridized with O 
2^ states. 

The LSDA calculations for BiCoOa suggest it is an insulator with an energy gap of 0.72 
eV. We obtained HS state of d^ electrons with magnetic moment of 2.41//^ which is smaller 
than 4 /is (expected for HS state) due to covalency effects. 

With the LSDA+U method BiNiOs is insulator with an energy gap of 1.23 eV. This value 
is twice as large as the experimental result from Reference 14^ but in good agreement with 
the present experimental estimation (see Table IVT) . The values of the calculated magnetic 
moments given in Table |V] are in agreement with the experimental one. From the density 
of states it can be seen that the top of the valence band and the bottom of the conduction 
band are formed by Bil 6s and Bi2 6s states hybridized with O 2p states. The occupancies 
calculated for the d states of the four non-equivalent Ni atoms in the LSDA+U calculation 
are 8.2, 8.36, 8.23, and 8.49 indicating that Ni has a 2+ valence which is in agreement with 
calculations performed in the literature^. 



B. Experimental results and discussion 

The O XES and O Is XAS measurements of BiMeOs (Me = Sc, Cr, Mn, Fe, Co, Ni), 
which probe the occupied and vacant O 2p states, respectively, are presented in Figure [3l 
The fine structure and energy distribution of the O Kq, XES matches the O 2^ occupied 
DOS obtained from LSDA calculations (see Figure [2]). The band gap for these materials 
was estimated using the peaks in the second derivative, this method has been shown to work 
well with O Kq, XES and Is XAS^^. A small amount of hybridization between the Bi 6s- 
and O ^j»-states is visible at about 519 eV in all compounds, as suggested by the calculated 
O 2p-stdXes. 

The estimated energy gaps from the experimental spectra and the calculations are not 
exactly the same. This is expected, since LSDA calculations are known to provide a un- 
derestimation of the energy gap, and the LSDA+U calculations used a typical value for U 
rather than a material-specific value. We have found that the calculated energy gaps are 
mostly within 0.5 eV of the experimental energy gaps and, more importantly, the trend in 
the calculated energy gaps is essentially the same as that in the experimental energy gaps 
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for these materials. The energy gaps are shown in Figure [H note the similar shape in the 
experimental and calculated curves. 
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FIG. 2: Calculated LSDA (LSDA-hU for BiMnOs and BiNiOs) partial density of states for BiMeOa 
(Me = Sc, Cr, Mn, Fe, Co, Ni) compounds. The DOS have been convoluted with the Fermi step 
at T = 300 K and broadened by 0.3 eV using a Lorentzian function to mimic the experimental 
resolution. The Fermi level corresponds to zero. (Colour in on-line version.) 

The comparison of the oxygen and iron X-ray spectra of BiFeOa with the calculated O 
2p and Fe 3d DOSes is shown in Figure O and Figure El respectively (solid line). The Fermi 
level position for the Fe Lg XES was determined using the XPS Fe 2p3/2 binding energy 
(710 eV)^^ In the absence of XPS O Is data the O K„ XES is compared with O 2p DOS by 
alignment of Bi 6s related subbands. The LSDA calculation reproduces the energy position 
of the Fe L3 XES (which probes occupied Fe 3d states) with respect to the Fermi level. 
On the other hand, the available LSDA-I-U calculations^!^ show a low-energy shift of Fe 3d 
band which contradicts the experimental Fe L3 XES. The spectral weight of Fe 3d- and O 
2p states is redistributed with increasing U value within LSDA-I-U calculations, as shown in 
Figures [5] and [61 As the U value rises the Fe 3d states are shifted to lower energies (Figure [6]). 
At the same time the O 2p states (Figure [5D are shifted towards Fermi level. Taking this 
trend into account both the occupied O 2p- and Fe 3d states obtained with U = eV are 
in better agreement with experimental spectra than the one calculated with a nonzero U. 
The Fe Lg^^ XES spectrum has some intensity at and just below the Fermi level while the 
O Kq, XES spectrum indicates a clear energy gap (see the lower panel of Figure ED in this 
energy region. In a recently published paper- it is concluded that LSDA-I-U calculations ( U 
= 2 eV) are a better match to experimental Fe L3 XES of BiFeOs than a LSDA (where 
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FIG. 3: Measured oxygen Kq, X-ray emission spectra (XES) and O Is X-ray absorption spectra 
(XAS) of BiMeOs compounds. The XAS is measured in total fluorescent yield (TFY) mode. The 
energy gap is estimated using the peaks in the second derivative, plotted below the corresponding 
spectra. The Bi 6s — O 2p hybridization peak at about 519 eV has been magnified by a factor of 
5. 




— • — Calculation 
— •— Experiment 



Sc Cr Mn Fe Co Ni 
Compound: BiMeO^ 

FIG. 4: The experimental (gray dots) and calculated (black dots) energy gaps for the BeMeOa 
compounds. 



U = eV) calculation. However, the comparison procedures used in matching calculations 
to experimental spectra are not described in Reference g]. Based on the comparison from 



Reference |6| one can estimate an XPS Fe 2p3/2 binding energy of 713 eV which contradicts 



15 




FIG. 5: The comparison between the O Kq XES and O Is XAS (lower panel) and the LSDA+U 
O 2p DOS (upper panel) of BiFeOa. In the upper panel, the calculated LSDA+U O 2p DOS 
for several different values of U are shown. The O 2p DOSes were broadened by 0.3 eV using a 
Lorentzian function to mimic the experimental resolution. The Fermi level corresponds to zero. 

the available experimental data (710 eV)^. 

As we have shown in Reference 5^, the O 2p DOS for binary oxides of non-transition 
elements occupy a whole energy range of excited d states of cations, involved in formation of 
the chemical bonding in solids. In BiMeOs compounds, as seen in Figure [2], a correspondence 
of energies in the main maxima of O 2p and Bi 6d DOS occurs. The main maxima of the 
calculated Me 3d DOS are also close to these energies. The non-bonding localized 3d states 
weakly mix with O 2p states near Fermi level. 

All O Kq, XES consist of low-energy subband located at ~2.5 eV with respect to the 
main maximum. As it was shown in Reference |56i, the low-energy subband of O XES of 
all binary oxides is directly connected with the s valence states of the cation. The energy 
separation of this subband from the main maximum depends in particular on the period 
number of Me. For ZnO it is ~ 2.4 eV, which coincides with the separation of subband from 
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FIG. 6: The comparison between the Fe Xg^j XES (lower panel) and the LSDA+U Fe 3d DOSes 
(upper panel) of BiFeOs. In the upper panel, the calculated LSDA+U Fe L2^3 DOS for several 
different values of U are shown. The Fe 3d DOS has been convoluted with the Fermi step at T 
= 300 K and broadend by 0.3 eV using a Lorentzian to mimic the experimental resolution. The 
Fermi level corresponds to zero. 

the main band of O K„ XES for BiMeOs compounds. The Bi 6s states are located ~ 5 eV 
below the Me 4s states, this is why they hybridize only weakly with the O 2p states (see 
Figure ED. 

The near edge fine structure of the O Is XAS (see Figure [3]) consists of three peaks 
(labeled as 1, 2 and 3) separated by ~ 1.5 - 2 eV for all compounds. There is a similar 
peak structure in the calculated unoccupied O 2p DOS (Figure [2]) in spite of the fact 
that the band gap (E^*^'^') is not reproduced in LSDA calculations. Indeed, the lack of 
the appropriate energy gap does not effect the correct reproduction of the energy distance 
between the main maximum of the O XES and the first peak of the O Is XAS (~ 5 
eV) in the LSDA+U calculations. This further demonstrates the well-known fact that O 2p 
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states are less influenced by correlation effects than d states—. 

The O Is XAS spectra of BiNiOa and BiCoOs show additional peaks near the Fermi 
level similar to those which were observed in the O Is XAS of LLNii_a;0^, LiCoO^^ and 
LaNiO^^. These additional peaks were attributed in References l57l]59l to the formation of O 
2p hole electronic states because of the energy disadvantage associated with a Ni^+ charge 
state. The Co^^ charge state is more stable, hence the intensity of the hole-like peak in 
BiCoOs is much lower than that in BiNiOa. Indeed, the "hole effect" in BiNiOs is even 
larger than that in Lij-Nii.^O^^. Among the materials studied in Reference 1571 the intensity 



of the O 2p hole peak was found to be largest for Lio.5Nio.5O where the hole peak had ~70% 
of the intensity of the main maximum of the O Is XAS. Similar intensity ratios have been 
reported for LaNiO^. In BiNiOs, however, the relative intensity of the O 2p hole peak is 
about 100% of the intensity of the main O 2p band. The O 2p holes are created in order 
to maintain the neutrality of the compound with the Ni^+ valence state. They are clearly 
observed in the O Is XAS for Ni oxides LaNiOs and BiNiOa. In the LSDA+U calculation 
the gap arises betwen the states predominantly formed by Bi 6s states strongly hybridized 
with O 2p states. The Ni atoms have a 2+ valency in agreement with experimental findings. 

Now let us discuss the origin of ferroelectricity in BiMnOs and BiFeOs. From Figure [2] 
it is clear that the Bi 6s states in all compounds are not hybridized with the Bi 6p states. 



As was supposed in Reference 



60, these Bi 6s lone pairs are slightly hybridized with O 



states, and can induce dipole moment. However, we could not find any difference in degree 
of Bi 6s-0 2p hybridization for BiMnOa and BiFeOs compared to that of other BiMeOs 
compounds where no ferroelectricity is observed. The stereochemical activity of the Bi lone 
pair induces the electric moment but if the symmetry of the system includes inversion the net 
polarization turns to zero. The space group for the monoclinic structure of BiMnOs contains 
inversion while the rhombohedral BiFeOs has R3c symmetry without inversion. Therefore 
the origin of ferroelectricity in these two compounds seems to be different. In BiMnOs the 
emergence of electric moment is connected with inversion symmetry breaking by the AFM 
order induced by particular orbital order below T = 474K— Another scenario for the 
onset of electric moment is possible in BiFeOs: an extensive investigation of ferroelectricity 
in BiFeOs in the framework of density functional theory was reported in Reference 
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It 



was mentioned that the driving force of the ferroelectric distortion is the Bi lone pair. The 
calculated values of polarization agree well with experimental results^. 
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Among the BiMeOs series there are three compounds with no inversion symmetry: 
BiMnOs (monochnic structure, but inversion is broken by AFM order), BiFeOs (rhom- 
bohedral structure) and BiCoOs (tetragonal structure). The ferroelectric properties of the 
BiMnOa and BiFeOs are confirmed experimentally while a ferroelectric hysteresis loop has 
not yet been observed for BiCoOs. One possible reason could be a low resistivity to the ap- 



plied electric field as was pointed out in Reference |4^. The crystal structures of the BiScOs 



{C2/c), BiCrOs {C2/c), and BiNiOs (PI) compounds include the inversion operation. 



V. SUMMARY 



In conclusion, the electronic structure of the BiMeOs (Me = Sc, Cr, Mn, Fe, Co, Ni) 
series was studied by soft X-ray emission and absorption spectroscopy. Experimental spectra 
were found to be in good agreement with spin-polarized electronic structure calculations. 
The presence of holes in the O Is XAS spectrum of BiNiOs was attributed to the 2+ valency 
of Ni. For all BiMeOa {Me = Sc, Cr, Mn, Fe, Co, Ni) compounds the band gap values were 
estimated from O XES and O Is XAS spectra. In case of BiFeOa this estimation results 
in the band gap value ~ 0.9 eV. The X-ray spectra for multiferroic and nonmultiferroic 
compounds in the series do not reveal any difference in degree of Bi 6s-0 2p hybridization. 
It was concluded that for nonzero electric polarization the stereochemical activity of the Bi 
lone pair should be accompanied with the inversion symmetry breaking. This condition is 
satisfied in case of BiMnOs, BiFeOs and BiCoOs but the electric polarization arises due to 
the different reasons. In BiMnOs noncolinear magnetic order with AFM component breaks 
the inversion symmetry and allow net polarization. In BiFeOs the rhombohedral crystal 
structure does not contain inversion and the crucial role is played by Bi lone pairs. BiCoOs 
is supposed to be ferroelectric but a ferroelectric hysteresis loop has not been observed as 
of now so apparently BiCoOs should be considered as a pyroelectric material. 
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